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1. Introduction 

The experimental evidence for discontinuities in the sequence of atomic masses, a— 
and j3- decay systematics and binding energies of nuclei suggests the existence of a set 
of magic proton and neutron numbers, which can be described successfully by single 
particle shell models with a heuristic spin-orbit term [T], [2J. The most prominent 
representative is the phenomenological Nilsson model j3] with an anisotropic oscillator 
potential: 



Although these models are flexible enough to reproduce the experimental results, they 
lack a deeper theoretical justification, which becomes obvious, when extrapolating the 
parameters k, /i, which determine the strength of the spin orbit and I term to the 
region of superheavy elements [5] . 

Hence it seems tempting to describe the experimental data with alternative 
methods. Typical examples are relativistic mean field theories [B], [7], where nucleons 
are described by the Dirac-equation and the interaction is mediated by mesons. 
Although a spin orbit force is obsolete in these models, different parametrizations 
predict different shell closures [8], [9]. Therefore the problem of a theoretical foundation 
of magic numbers remains an open question since Elsasser |10j raised the problem 75 
years ago. 

A fundamental understanding of magic numbers for protons and neutrons may be 
achieved if the underlying corresponding symmetry of the nuclear many body system 
is determined. Therefore a group theoretical approach seems appropriate. 




(1) 
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Group theoretical methods have been successfully applied to problems in nuclear 
physics for decades. Elliott [IT] has demonstrated, that an average nuclear potential 
given by a three dimensional harmonic oscillator corresponds to a SU(3) symmetry. 
Low lying collective states have been successfully described within the IBM-model 
[T2], which contains as one limit the five dimensional harmonic oscillator, which is 
directly related to the Bohr-Mottelson Hamiltonian. 

In this paper we will determine the symmetry group, which generates a single 
particle spectrum similar to (JT|), but includes the magic numbers right from the 
beginning. 

Our approach is based on group theoretical methods developed within the 
framework of fractional calculus. 

The fractional calculus [13]- [16] provides a set of axioms and methods to extend 
the coordinate and corresponding derivative definitions in a reasonable way from 
integer order n to arbitrary order a: 

gn Qa 

(*"^>^^> < 2 > 

The definition of the fractional order derivative is not unique, several definitions e.g. 
the Feller, Fourier, Riemann, Caputo, Weyl, Riesz, Griinwald fractional derivative 
definitions coexist [T7]-[5!)]. A direct consequence of this diversity is the fact, that the 
solutions e.g. of a one dimensional wave equation differ significantly depending on the 
specific choice of a fractional derivative definition. 

Until now it has always been assumed, that the fractional derivative type for an 
extension of a fractional differential equation to multi-dimensional space should be 
chosen uniquely 

In contrast to this assumption, in this paper we will investigate properties of 
higher dimensional rotation groups with mixed Caputo and Riemann type definition 
of the fractional derivative. We will demonstrate, that a fundamental dynamic 
symmetry is established, which determines the magic numbers for protons and neutron 
respectively and furthermore describes the ground state properties of nuclei with 
reasonable accuracy. 

2. Notation 

We will investigate the spectrum of multi dimensional fractional rotation groups for 
two different definitions of the fractional derivative, namely the Riemann- and Caputo 
fractional derivative. Both types are strongly related. 

Starting with the definition of the fractional Riemann integral 

um)(x) = ~ r d £ ( x - o a - x /(o x > o 



+J M ' T(a) J n 

ul a f(x)={ \' J0 (3) 

Ul ~ f){x)= T(ajJ *^- a; ) a ~V(0 *<0 

where T(z) denotes the Euler T-function, the fractional Riemann derivative is defined 
as the product of a fractional integration followed by an ordinary differentiation: 

B a — — / 1_a (A) 
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cBS = ^ 1 ~ a iz (6) 



r(l - a) dx J n 

The Caputo definition of a fractional derivative follows an inverted sequence of 
operations Q. An ordinary differentiation is followed by a fractional integration 

d_ 
dx 

which results in: 

^/(*)=<j \ ° „ (7) 

Applied to a function set f(x) — x na using the Ricmann fractional derivative definition 
(|5| we obtain: 

fl« t «q _ T(l + no) (n-l)a ro\ 
R ° xX "r(l + (n-l)a) X W 

= K [n]x< n - r > a (9) 

where we have introduced the abbreviation R [n\. 

For the Caputo definition of the fractional derivative it follows for the same 
function set: 

f r ( 1 + nQ ) x (n-l)« „ > 

c ^x ,lQ = I r(l + (n-l)a) " (10) 

1 n = 

= c [n] a; (n - 1)Q (11) 

where we have introduced the abbreviation c [n]. 

Both derivative definitions only differ in the case n = 0: 

cN=rN-(S„ok[0] (12) 

= rN - S n0 r _ (13) 



where <5 mn denotes the Kronecker-(5. We will rewrite equations ([9| and (111 
simultaneously, introducing the short hand notation 

R , c d^x na = R , a [n]x^ a (14) 

We now introduce the fractional angular momentum operators or generators of 
infinitesimal rotations in the i,j plane on the iV-dimcnsional Euclidean space: 

R , c Lij(a) = H>c ih(x?d? - x?8?) (15) 

The commutation relations of the fractional angular momentum operators are 
isomorph to the fractional extension of the rotational group SO(N) 

r,c [Lij(a), L k i(a)} = RjC ihf i]k i mn L mn {a) (16) 

i,j,k,l,m,n = 1,2,.., N 




Figure 1. Spectrum of the Casimir operator L 2 (L, a) from \l9\ as a function of 
the fractional derivative coefficient a. Only the L = state differs for Riemann 
and Caputo derivative. 



with structure coefficients R , c fijki mn ■ Their explicit form depends on the function set 
the fractional angular momentum operators act on and on the fractional derivative 
type used. 

The Casimir-operators of the fractional rotation group R SO a (3) based on the 
Riemann fractional derivative definition have been derived in |26] and for c SO a (3) 
based on the Caputo fractional derivative definition are given in |27j . We summarize 
the major results: 

According to the group chain 

R , c SO a (3) D RC SO a (2) (17) 
there are two Casimir-operators Aj, namely A2 = L z (a) = £12(0:) and A3 = L 2 (a) = 
L12 ( a ) + -^13 (ct) + L\ z (ot). We introduce the two quantum numbers L and M, which 
completely determine the eigenfunctions \LM >. It follows 

R , c L z {a)\LM > = RjC fisign(M) [\M\] \LM > (18) 

M = —L, —L + 1, ±0, L 

R _ c L 2 (a)\LM> = R , c h?[L][L + l]\LM > (19) 

L = 0,l,2,... 

where \M\ denotes the absolute value of M. In addition, on the set of eigenfunctions 
\LM >, the parity operator II is diagonal and has the eigenvalues 

U\LM >=(-l) L \LM > (20) 
In figure [l] the eigenvalues of the Casimir-operator L 2 are shown as a function of a. 
Only in the case L = the spectra differ for the Riemann- and Caputo derivative. 
While for the Caputo derivative 

c £ 2 (a)|00 >= (21) 

because c [0] = 0, using the Riemann derivative for a ^ 1 there is a nonvanishing 
contribution 

R L 2 (a)\00 >= R ft 2 [0][l]|00 >= h 2 l^ + a \ \00 > (22) 

r(l-a) 
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We now use group theoretical methods to construct higher dimensional representations 
of the fractional rotation groups R:C SO a (S). 

As an example of physical relevance we introduce the group CRE G with the 
following chain of subalgebras: 

CRR G D c SO a (3) D R SO Q (3) D R SO a (3) (23) 

The Hamiltonian H can now be written in terms of the Casimir operators of 
the algebras appearing in the chain and can be analytically diagonalized in the 
corresponding basis. The Hamiltonian is: 

H=^ c Ll(a) + ^Lj(a) + ^ n Ll(a) (24) 

with the free parameters W\,u) 2 , ^3 and the basis is \LiM\L 2 M 2 L 3 M 3 >. Furthermore, 
we impose the following symmetries: 

First, the wave functions should be invariant under parity transformations, which 
according to ( 20 1 leads to the conditions 

L x = 2ni L 2 = 2n 2 L 3 = 2n 3 , m, n 2 , n 3 = 0, 1, 2, 3, ... (25) 

second, we require 

c L Zl (a)\LiMiL 2 M 2 L 3 M 3 > = c + H[Li]\LiMiL 2 M 2 L 3 M 3 > (26) 

R L Z2 (a)\L 1 M 1 L 2 M 2 L 3 M 3 > = R + H[L 2 ]\L 1 MiL 2 M 2 L 3 M 3 > (27) 

H L Z3 {a)\LiMiL 2 M 2 L 3 M 3 > = R + ft[L 3 ] |£iMi L 2 M 2 L 3 M 3 > (28) 

which reduces the multiplicity of a given \2niMi2n 2 M 2 2n 3 M 3 > set to 1. 

With these conditions, the eigenvalues of the Hamiltonian (24 1 are given as 

E(a) = hu lc [2ni][2ni + 1] + 

tkJ 2H [2n 2 ] [2n 2 + 1] + huo 3R {2n 3 ] [2n 3 + 1] (29) 

A r(l + (2n t + l)a) r(l + a) 

= lj^r(l + (2n i -l)a) ~ 6n > otUJl r0^a) (30) 

ni,n 2 ,n 3 = 0,1,2,.. 

on a basis \2ni2ni2n 2 2n 2 2n 3 2n 3 >. 

This is the major result of our derivation. We call this model the Caputo- 
Riemann-Riemann symmetric rotor. What makes this model remarkable is its 
behaviour near a — 1/2. 

On the left of figure [2] we have plotted the energy levels in the vicinity of a ps 1/2 
for the case 

^1=^2 = ^3= UJq (31) 

which we denote as the spherical case. 

For the idealized case a = 1/2, using the relation T(l + z) = zT(z) the level 
spectrum ( 30 1 is simply given by: 

3 1 

E(a = 1/2) = &Jo(ni + n 2 + n 3 + - - -S ni0 ) (32) 

For ni 7^ this is the well known spectrum of the 3-dimensional harmonic oscillator. 
Assuming a twofold spin degeneracy of the energy levels, we introduce the quantum 
number TV as 

N = m + n 2 + n 3 (33) 



Fractional dynamic symmetries and the ground state properties of nuclei 6 




Figure 2. On the left the energy spectrum E{a) from ( | 30 [ > for the spherical 
case ( |31| in units of Hluq for the Caputo-Riemann-Riemann symmetric rotor is 
presented. The right diagram shows the neutron energy levels for the spherical 
nucleus \ \ \X calculated within the framework of the asymmetric two center shell 
model (ATCSM) 1281 near the ground state as a function of increasing strength 
of the spin-orbit term (kokIs) increasing from 50% to 100% of the recommended 
ft value, while the fil 2 value is kept constant. The transition from magic numbers 
of the standard 3-dimensional harmonic oscillator levels j34} to the shifted set of 
magic numbers l |37| l is pointed out with thick lines. Left and right figure therefore 
show a similar behaviour for the energy levels. 



Consequently we obtain a first set n ma gj c i of magic numbers ™ ma gj c 

^magicl = ^ + 1 )( iV + 2 )(^+ 3 ) W= 1,2,3,... (34) 

= 8,20,40,70,112,168,240,... (35) 

which correspond to the standard 3-dimensional harmonic oscillator at energies 

E = hv (N + 3/2) (36) 

In addition, for n\ = 0, which corresponds to the 1 00 In^ln^ 2n^2ns > states, we 
obtain a second set ^ ma gj c 2 °f magic numbers 

"magic 2= ^+l)(5 + (iV+l) 2 ) N = 0, 1, 2, 3, ... (37) 

= 2, 6, 14, 28, 50, 82, 126, 184, 258, ... (38) 

at energies 

E = huj (N + l) (39) 

which is shifted by the amount — -^Hu>o compared to the standard 3-dimensional 
harmonic oscillator values. 
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From figure p| it follows, that for a < 1/2 the second set « ma gi c 2 °f energy levels 
falls off more rapidly than the levels of set ^ ma gj c \ ■ As a consequence for decreasing 
a the magic numbers ^ ma gi c \ die out successively. On the other hand, for a > 1/2 
the same effect causes the magic numbers " ma gi c \ to survive. 

We want to emphasize, that the described behaviour for the energy levels in 
the region a < 1/2 may be directly compared to the influence of a ^s-term in 
phcnomenological shell models. As an example, on the right hand side of figure [2] 
a sequence of neutron levels for the superheavy element ffl-X" calculated with the 
asymmetric two center shell model (ATCSM) [25] with increasing strength of the Is- 
term from 50% to 100% is plotted. It shows, that the n = 168 gap breaks down at 
about 70% and the n— 112 gap at about 90% of the recommended k- value for the Is- 
term. This corresponds to an a *=s 0.46 value, since in the Caputo-Riemann-Riemann 
symmetric rotor the n — 168 gap breaks down at a = 0.466, the n — 112 gap at 
a = 0.460 and the n — 70 gap vanishes at a = 0.453. 

We conclude, that the Caputo-Riemann-Riemann symmetric rotor predicts a well 
defined set of magic numbers. This set is a direct consequence of the underlying 
dynamic symmetries of the three fractional rotation groups involved. It is indeed 
remarkable, that the same set of magic numbers is realized in nature as magic proton 
and neutron numbers. 

In the next section we will demonstrate, that the proposed analytical model is an 
appropriate tool to describe the ground state properties of nuclei. 



4. Ground state properties of nuclei 



We will use the Caputo-Riemann-Riemann symmetric rotor ( 30 1 as a dynamic shell 
model for a description of the microscopic part of the total energy E^ of the nucleus. 

-^tot = ^macroscopic + ^microscopic (40) 
= ^"macroscopic + + ^ (41) 

where SU and SP denote the shell- and pairing energy contributions. 

For the macroscopic contribution we use the finite range liquid drop model 
(FRLDM) proposed by Moller [31] using the original parameters, except the value 
for the constant energy contribution ao- 

As the primary deformation parameter we use the ellipsoidal deformation Q: 

Q= b - = ^ = ^ (42) 
a u>i u>2 

where a, b are the semi-axes of a rotational symmetric ellipsoid. Consequently a value 
Q < 1 describes prolate and a value of Q > 1 descibes oblate shapes. In order to 
relate the ellipsiodal deformation Q to the quadrupole deformation 62 used by Moller, 
we define: 

0= 1 - 1.43085e 2 + 0.707669e2 (43) 

Furthermore we extend the original FRLDM-model introducing an additional 
curvature energy term Vr{Q), which describes the interaction of the nucleus with 
the collective curved coordinate space [32] : 

V R (Q) = ~a R B R A- 5 ^ (44) 
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Figure 3. As a test of the plateau condition dU/d-f = for the Strutinsky shell 
correction method, the total shell correction energy SU = SU p + SUn is plotted 
as a function of 7 for different nuclei. 



where A is the nucleon number, an is the curvature parameter given in [MeV] and the 
relative curvature energy B R (Q) given as: 

R mi-qnW 199 - 2881og(2) \ 2 

MW ^ ^ (2 + Q2)(266 _ 67Q2 + 96(g2 _ 4)log(2)) j 

which is normalized relative to a sphere B R (Q = 1) = 1. 
Therefore the total energy may be splitted into 

E tot = £W + E mic (46) 

where 

E mac (a , Ofl) = FRLDM(oo, 0=1) + V R (Q = 1, a R ) (47) 
^ m ic(ao, a R , Q) = +5U + 6P + FRLDM(a , Q) + V R (Q, a R ) 

- (FRLDM(a , Q = 1) + Vr(Q = 1, a«)) (48) 

with two free parameters ao,a R , which will be used for a least square fit with the 
experimental data. 

For calculation of the shell corrections we use the Strutinsky method [29 ,[30 . 
Since we expect that the shell corrections are the dominant contribution to the 
microscopic energy, for a first comparison with experimental data we will neglect 
the pairing energy term. 

In order to calculate the shell corrections, we introduce the following parameters: 



Hu)q 


= 38 A- i [MeV] 




(49) 








(50) 


L02 


= U3\ 




(51) 


LJ 3 






(52) 






' 0.46 + 0.000220 Z 


Z > 50 




OLZ 


-\ 


0.2469 + 0.00448 Z 


28 < Z < 50 


(53) 






k 0.2793 + 0.00332 Z 


Z < 28 








' 0.41 + 0.000200 N 


N > 50 




a N 


■\ 


0.3118 + 0.00216 N 


28 < N < 50 


(54) 






^ 0.2793 + 0.00332 N 


N < 28 
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Figure 4. Comparison of calculated shell corrections SU from the Caputo- 
Riemann-Riemann symmetric rotor l ]30| with the parameter set l |49[ l-| [56] l (thick 
line) with the tabulated E m ; c from Moller |31| (dashed line). Upper two rows 
show values for a given Z as a function of N, lower two rows for a given N as a 
function of Z. Bars indicate the experimentally known region. The original 62 
values from 1311 are used, which is the main source of error. 



7 
m 

«o 
an 



= 1.1 Hluq 

= 4 

= 2.409[MeV] 
= 15.0[MeV] 



(55) 
(56) 
(57) 
(58) 



Input parameters are the number of protons Z, number of neutrons N, the nucleon 
number A — N + Z, and the ground state quadrupole deformation ti. 

The values obtained include the frequencies (50 1, (51 1, (52 1, which result from a 
least square fit and quadratic approximation of equipotcntial surfaces, the fractional 



derivative coefficients for protons ( 53 1 and neutrons ( 54 ) which determine the level 



spectrum for protons and neutrons for the proton and neutron part of the shell 
correction energy respectively from a fit of the set of nuclids ll^i, 5o°Sn 
and from the requirement, that the neutron shell correction for £n Sn should amount 



132c 
50 



208T3U 

82 rD 



about — 5.1[MeV], (551 from the plateau condition dU/dj = (see figure [3]) and (561 
the order of included Hermite polynomials for the Strutinski shell correction method. 
Finally Itioq, a , an from a fit of the experimental mass excess given in |38j . 

We compare our results with for the microscopic energy contribution -E m } c with 
data from Moller et. al. [3T] and use their tabulated €2 values. They have not only 
listed data for experimental masses but also predictions for regions, not yet confirmed 
by experiment. 



Fractional dynamic symmetries and the ground state properties of nuclei 



10 



0.2 



0.0' 



-0.2 
5 



5 
1- 



-10 
5 



,_, 
L±J 

-10 



> 5 
I 

UJ -5 



1 1 1 1 1 [ 1 . 1 1 [ 1 1 1 [ 1 1 

■...-,» . ! ! 1 ■■_ .i 1 .■■■■".ii:- „. j ■■- ■—■■ llB i-. „ l' 1 c 2 
_ ;j I./ V ■■i" 5 " 1 "- !,,i:S " 

■ ^■■^H W crrG 

- , 1 ■ 1 ■ 1 . 1 , 1 T 1 i 1 i T 


1 1 1 1 


1 ■ 1 1 1 1 1 1 1 1 1 _ 

^^^^^^^^^^^^^^^^ ^ m ' c th 




: M 


-+ I 1 1 1 1 1 1 1 1 h 

^^L^ Emic 


1 1 


-j^J 1 1 1 1 1 1 1 1 h 



20 



40 



60 



80 

N 



100 



120 



140 



160 



Figure 5. Comparison of calculated i? m j c from the Caputo-Riemann-Riemann 
symmetric rotor |30| with the parameter set |49|-|56]l, minimized with respect 
to €2 with the experimental masses from Audi ]38 as a function of N. From top 
to bottom the minimized £2 values, theoretical £? m ; c , experimental microscopic 
contribution from the difference of experimental mass excess and and macroscopic 
FRLDM energy and error in [MeV] are plotted. 



In figure [4] we compare the calculated 8U values with the tabulated -E m j c , which 
is justified for almost spherical shapes (e 2 ~ 0). The results agree very well within the 
expected errors (which are estimated ss 2 [MeV] for the pairing energy and 0.5 [MeV] 
for E m ^ c ), especially in the region of experimentally known nuclei. 

A remarkable difference between the calculated shell correction and tabulated 
-E m } c from Moller occurs for superheavy elements (N — 184, last picture in figure 
[4|. While phenomenological shell models predict a pronounced minimum in the shell 
correction energy for Z — 114 [33]- [37J the situation is quite different for the rotor 
model, where two magic shell closures at Z = 112 and Z — 126 are given, but the 
Z = 112 shell closure is not strong enough to produce a local minimum in the shell 
correction energy plot as a function of Z. Instead, between Z = 112 and Z — 126, 
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there emerges a slightly falling energy plateau, which makes the full region promising 
candidates for stable, long-lived superheavy elements. 

While this result contradicts predictions made with phenomenological shell 
models, it supports recent results obtained with relativistic mean field models [7J, 
which predict a similar behaviour in the region of super heavy elements as the proposed 
rotor model. 

In figure [5] we have covered the complete region of available experimental data 
for nuclids and compare the calculated theoretical microscopic energy contribution 
minimized with respect to the deformation with the experimental values. The influence 
of shell closures is very clear. The rms-error is about 2.4[MeV]. The maximum 
deviation occures between closed magic shells. Therefore in the next section we will 
introduce a generalization of the proposed fractional rotor model, which not only 
determines the magic numbers accurately but in addition determines the fine structure 
of the single particle spectrum correctly. 



5. Fine structure of the single particle spectrum - the extended 
Caputo-Riemann-Riemann symmetric rotor 

In the previous section we have demonstrated, that the Caputo-Riemann-Riemann 
symmetric rotor correctly determines the magic numbers in the single particle spectra 
for neutrons and protons. However, there remains a significant difference between 
calculated and experimental ground state masses for nuclei with nucleon numbers 
far from magic shell closures. This is a strong indication for the fact, that the fine 
structure of the single particle levels is not yet correctly reproduced. 

We therefore propose the following generalization of the Caputo-Riemann- 
Riemann symmetric rotor group: 

C3C2R3R3 

G D c SO a (3) D c SO a (2) D R SO a (3) D R SO a (3) (59) 



with the Casimir operators ( 18 I and (h 9k it follows for the Hamiltonian H 



H = ^ L?(a) + Blo 0c L Zi (a) + ^ R L 2 2 (a) + —^(a) (60) 

with the free parameters u>\, u>2>i^3, B, where B may be called fractional magnetic field 
strength in units [fiwo]- 



Imposing the same symmetries (25), (26 1 as in the case of the symmetric Caputo- 



Riemann-Riemann rotor, the eigenvalues of the Hamiltonian (60) are given as 
E(a) = /kJi c [2ni][2ni + 1] + Btkj c[2ni] + 

/icj 2R [2n 2 ] [2n 2 + 1] + hu} 3R [2n 3 ] [2n 3 + 1] (61) 

A r(l + {2m + l)a) . , T(l + a) 
•f— ^ -1(1 + (2ftj — lja) I (1 — a) 

r(i + (gniH i 



T(l + (2m - l)o) 1U u T(l-a) 

Hi, 712,713 = 0,1,2, .. (62) 

on a basis |2ni2ni2n22ri22n32n3 >. 

We call this model the extended Caputo-Riemann-Riemann symmetric rotor. 
The additional a L Zl (a) term yields a level splitting of the harmonic oscillator set of 



magic numbers ^ ma gj c \ ( 34 ) , while the multiplicity of the «- ma gj c 2 se ^ ( 37 1 remains 




Figure 6. For a = 1/2, on the left side the level spectrum for the extended 
Caputo-Riemann-Riemann symmetric rotor \62\ is plotted as a function of the 
fractional magnetic field strength B. The levels are labeled according to the 
corresponding [Nlj z ] Nilsson scheme and the multiplicity is given. For a = 0.48 
the resulting level sequence near N « 126 is plotted on the right. At B fa 0.25 
the resulting spectrum coincides with the corresponding spherical Nilsson level 
spectrum. Brackets indicate the proposed appropriate combinations of rotor 
levels. 



unchanged, since this set is characterized by n\ = 0. This is exactly the behaviour 
needed to describe the experimentally observed fine structure, as can be deduced from 
the right hand side of figure [2| 

In order to clearly demonstrate the influence of the additional term, we first 
investigate the level spectrum for the spherical (31 1 and idealized case a = 1/2. 



The level spectrum (62 1 simply results as: 



1 



E(a = 1/2) = huj (ni + n 2 + n 3 + ^ - ^ n ,o) 
+ BHuj ( 



nil 
T(l/2 4 



1 



n 2 + n 3 4 
2 ni ni! 



r(l/2) 



(63) 



2 

^nio) 



(64) 



= Huj (ni + 

BhLUn 

-I -I 

^ l (2ni-l)!! 

where !! denotes the double factorial. 

On the left side of figure [6]this spectrum is plotted in units [Hujq]. Single levels are 
labeled according to the Nilsson-schcme and multiplicities are given in brackets. For 
small fractional field strength B the resulting spectrum exactly follows the schematic 
level diagram of a phcnomcnological shell model with spin-orbit term, as demonstrated 
e.g. by Goeppert-Mayer [T]. 

A small deviation from the ideal a = 1/2 value reproduces the experimental 
spectra accurately: For a = 0.48 the resulting level spectrum is given on the right hand 
side of figure[6] Obviously there is an interference of two effects: First, for a =/= 1 /2 now 
the degenerated levels of both magic sets split up and second the fractional magnetic 
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Figure 7. Comparison of cal culated S m ; c from the ex ten d ed C aputo-Riemann- 
Riemann symmetric rotor | |62[ l with the parameter set (|65|-l |68| l, minimized with 
respect to E2 with the experimental masses from Audi [38] as a function of N. 
From top to bottom the minimized €2 values, theoretical and experimental masses 
and error in [MeV] are plotted. 



field B acts on the subset ^ ma gi c \- For B « 0.25 the spectrum may be directly 
compared with the spherical Nilsson level scheme, which is given for neutrons between 
82 < N < 126 as 2/|, 1/if, H^, 3p§, 2/f, 3p|, see e.g. results of [39], which 
corresponds to a sequence of sub-shells at 90, 100, 114, 118, 124, 126. This sequence is 
correctly reproduced with the extended Caputo-Riemann-Riemann symmetric rotor. 

With the parameter set, which is obtained by a fit with the experimental masses 
of Ca-, Sn- and Pb-isotopes 

huj = 28 A- ^ [MeV] (65) 

f 0.480 + 0.00022 Z Z > 50 
az = < 66 
I 0.324 + 0.00332 Z Z < 50 
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I 0.446 + 0.00022 N N > 29 
a N = I (67) 
1 0.356 + 0.00332 TV N < 29 y ' 

B = 0.27 (68) 

the experimental masses are reproduced with an rms-error of 1.7[MeV]. Results are 
given in figure [7] 

The deformation parameters, obtained by minimization of the total energy, are to 
a large extend consistent with values given in [3T] e.g. for 264 Hsi 8 we obtain e 2 = 0.22, 
which conforms with Moller's (e 2 = 0.2) and Rutz's results [5]. However, there occur 
discrepancies mostly for exotic nuclei. For example our calculations determine the 
nucleus 42 Si to be almost spherical, while Moller predicts a definitely oblate shape. 

Finally, defining a nucleus with e 2 > 0.05 as prolate and with e 2 < —0.05 as oblate 
the amount of prolate shapes is about 74% of all deformed nuclei. This is close to the 
value of 82% [40] . obtained with the Nilsson model using the standard parameters. 

Summarizing the results presented, the proposed extended Caputo-Riemann- 
Ricmann symmetric rotor describes the ground state properties of nuclei with 
reasonable accuracy. We have demonstrated, that the nuclear shell structure may 
indeed be successfully described on the basis of a dynamical symmetry model. 

The advantages of this model, compared to phcnomcnological shell and relativistic 
mean field models respectively are obvious: 

Magic numbers are predicted, they are not the result of a fit with a 
phcnomcnological Zs-term. There are no potential-terms or parametrized Skyrme- 
forces involved and finally, all results may be calculated analytically. 

The results obtained encourage further investigations in this field. The next 
steps should include the pairing energy term and parameters should be determined 
by a more sophisticated fit procedure. With these additional contributions the model 
will most probably describe nuclear properties with at least similar accuracy as the 
models currently used. 

6. Conclusion 

Based on the Riemann- and Caputo definition of the fractional derivative we used 
the fractional extensions of the standard rotation group SO(3) to construct a higher 
dimensional representation of a fractional rotation group with mixed derivative types. 

We obtained an extended symmetric rotor model, which predicts the sequence of 
magic proton and neutron numbers accurately. Furthermore we have shown, that the 
ground state properties of nuclei can be reproduced correctly within the framework of 
this model. 

Hence we have demonstrated, that a dynamic symmetry, generated by mixed 
fractional type rotation groups is indeed realized in nature. 
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